Measurement of the vortex mass in a superconducting film. 
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We have combined high resolution magneto-optical imaging with an ultra- fast heating/cooling 
technique to measure the movement of individual vortices in a superconducting film. The motion 
took place while the film was heated close to Tc, where pinning and viscous forces are relatively 
small. Under these conditions, vortices move due to the magnetic repulsion between them. We found 
that a finite vortex mass has to be included in the analysis in order to account for the experimental 
results. The extent of the motion is consistent with a vortex mass being 3 orders of magnitude 
smaller than the mass of all the electrons in the core. 



Vortices in a superconductor are localized, topologi- 
cally protected excitations carrying a quantized amount 
of magnetic flux. They are the "elementary particles" 
in several models of statistical mechanics(Kosterlitz- 
Thouless, XY model). One important parameter char- 
acterizing a particle is its mass. The size of the mass de- 
cides whether vortices are quantum or classical objects^. 
Quantum vortices in unconventional superconductors are 
predicted to support Majorana excitations^ and can be 
applied for quantum computing^. A quantized vortex 
has no intrinsic mass. An effective mass however can be 
associated with its motion. Due to its importance, the 
concept of the vortex mass was discussed extensively over 
the years but remains a controversial issue. One point of 
view is that mass plays no role in the dynamics since 
an inertial term in the equation of motion of a vortex is 
always negligible next to the viscous drag force. Predic- 
tions for the value of the mass in different limits (dirty 
and superclean) span 5 orders of magnitude^^— . Exper- 
imentally, the presence of an inertial term is difficult to 
detect since at low temperatures vortices in superconduc- 
tors are pinned, and if they move at all, their motion is 
dominated by viscosity. We are aware of only one such 
attempt, with inconclusive results^. It is quite clear that 
in order to check if an inertial term plays a significant role 
in vortex dynamics, one should realize conditions where 
the viscous drag force is not dominant. 

Both pinning^ and viscous drag forcesiS decrease 
strongly with temperature and vanish at the critical tem- 
perature Tc- In a thin superconducting film, the vortex- 
vortex interaction (Fi„t)is mainly due to the magnetic 
field of the vortices outside the filmii. Fi„t depends on 
temperature rather weakly. Consequently, a temperature 
interval exists near in which the viscous and pinning 
forces can be smaller than the vortex-vortex interaction, 
and do not dominate the dynamics. Our experiment was 
designed to measure vortex motion in this temperature 
interval. To obtain quantitative estimates, we use a clas- 
sical equation of motion of a vortex-S'^^. This equation 
is derived by integrating over the microscopic degrees of 
freedom, leaving only macroscopic forces: 



At„d— = F„t + Fi - VC/p - 771; (1) 



Here ^l^ is the vortex mass per unit length, d is 
the film thickness, Fi„t is the vortex-vortex interaction, 
Fi oc J X $0 is the Lorentz force due to interaction with 
currents J (<f>o is the flux quantum). Up is the pinning po- 
tential, and ryv is the viscous drag force. Measurements^^ 
show that in Nb, the Magnus forcoi^ is much smaller 
then the viscous force and can be neglected. In the ab- 
sence of external currents and fields, the Lorentz force 
results from currents associated with vortices trapped in 
the sample. Our observations were performed in a region 
near the center of the sample where the currents vanish 
and this force is small^^. 

The vortex-vortex interactionii Fi„t is repulsive for 
vortices of the same polarity. In the limit r ^ A (A is 
the penetration depth), Fi„t depends on the intervortex 
spacing r as Fint = ^0/ Hot^^^ ■ This relation holds even 
at O.OTTc, the maximal temperature where we use it. At 
this temperature A ~ 0.3/zm while the minimal r l/im. 

We determined the pinning and the viscous drag forces 
for our films from transport measurements using the 
method described by Klein et ali^. The temperature de- 
pendence of the various forces is shown in Fig. [TJ In the 
interval between T* w 0.95Tc and Tc F^t is larger than 
the viscous drag force. Our experiment was performed 
in this interval. 

The sample is a 200nm thick Niobium film with Tc 
of 8.8K deposited on a sapphire substrate. The film is 
patterned into squares of AOOfim x 400/iTO. The experi- 
ment is shielded by //-metal from external magnetic fields. 
Small fields (< mT) are applied using a solenoid inside 
the shield. The vortices in the film are imaged using 
high resolution magneto opticsi^. Using 10 s integration 
time, our system images large areas of 100 x lOO^m^ 
with a O.S^m resolution. Currently, our system shows 
the highest resolution achieved with magneto optics. 

Since the typical velocity of a vortejsi^ is few km/sec, 
it traverses the field of view in several nanoseconds. No 
technique can image this motion in real time. Our ap- 
proach is to take a pinned vortex array, release the pin- 
ning for about 1 nsec, allowing the vortices to move under 
the various forces, and then restore the pinning. Pinning 
can be turned off and on by rapidly heating and cool- 
ing the film, as illustrated in the inset of Fig. [TJ Dur- 
ing the short Infra-Red laser pulse, the temperature of 
the film is increased from the base temperature of 5.5K 
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FIG. 1. Temperature dependence of the forces acting on a vor- 
tex. The quantitative values are for a 200 nm Nb film. Solid 
black line is Fi„t for two vortices separated by Circles 
(red) show the maximal drag force taken from our transport 
measurements. Dashed line (red) is a linear fit. Blue dia- 
monds and dotted line represent the measured pinning force 
and the fit to Fp{T) = Fp(0)(l - T/Tcf. T" is the crossover 
temperature above which the pinning and the drag force are 
weaker than Fint and vortices are able to move. The inset 
shows the time dependence of the temperature of the film 
during a heating pulse. Motion of vortices takes place only 
during the short time interval where the temperature exceeds 
T*. 



by an amount proportional to the intensity of the light. 
The 1mm thick sapphire substrate is transparent at the 
wavelength of the laser. Hence, by back side illumination, 
only the film heats up, while the substrate remains near 
the base temperature. Heat escapes from the film via 
ballistic phonons crossing into the cold substrate. The 
substrate, which acts as a heat sink has a thermal mass 
about 1000 times larger than that of the film, and so its 
temperature does not increases significantly during the 
pulse. The thermalization time constant of the film is 
~ 6 • 10~^^sec, so essentially the film remains hot as long 
as the laser pulse is on, with its instantaneous temper- 
ature proportional to the light intensity. We confirmed 
this scenario by comparing the temperature profile mea- 
sured directly using a GeAu thin film bolometer. The 
laser pulse intensity profile was measured simultaneously 
using a fast photodiode. In our apparatus the length 
of the heating pulse is fixed, producing a temperature 
profile shown in the inset of Fig.l. The static positions 
of vortices are imaged before and after the motion took 
place. 

Our initial vortex configuration contains small aggre- 
gates only few vortices each, shown in Fig. [51 To prepare 
such aggregates, we first cool the film in a very small 
field, typically 30 ^T. The field is turned off at low tem- 
perature, leaving vortices trapped in the film. We then 
apply an inhomogeneous heating light pulse. The speckle 
pattern of the light creates an inhomogeneous illumina- 




FIG. 2. Typical Magneto-Optical image of vortex aggregates 
in the superconducting film. Each white spot is an individ- 
ual vortex. Inset: Light intensity distribution of a typical 
inhomogeneous laser heating pulse. Red spots represent high 
intensity regions, where the temperature during the pulse is 
above Tc- 



tion, strong enough to increase the temperature above 
Tc in some small regions. Some of the vortices escape 
into these normal regions and become trapped as classical 
fiux. Upon cooling, these regions become superconduct- 
ing again and the trapped fiux disintegrates into small 
aggregates of several vortices each. The typical distance 
between nearby vortices in an aggregate is w l^m. Typ- 
ically, the aggregates are separated by ~ 10/i-m, with the 
area in between largely free of vortices. 

We next apply a homogeneous heating pulse to the 
film containing the initial vortex array shown in Fig. [21 
To obtain uniform illumination (and heating), the light 
was passed trough a diffuser. The random intensity vari- 
ability across the sample translated to temperature units 
was less then O.OlTc. The displacement of vortices can 
be seen in the differential image shown in Fig. [3l gen- 
erated by subtracting the image before the pulse from 
the image after. If the maximal temperature during the 
pulse does not reach ^ 0.95Tc, we see no motion within 
our resolution. Once the maximal temperature during 
the pulse exceeded ^ O.QSTc, the aggregates begin to dis- 
perse (shown in Fig. [3^). This temperature is in good 
agreement with T* estimated using our transport mea- 
surements (see Fig. [J). At T = 0.97Tc, the typical dis- 
placement of the vortices is several /im (see Fig.[5l3). At 
T — O.OOTc, (Fig. [St) the displacement becomes compa- 
rable to the distance between aggregates. All the trapped 
vortices disappear at Tc. 

An important check involves the displacement of iso- 
lated vortices. At T„iax = O.OSTc isolated vortices do not 
move. At Tmax — 0.97Tc (Fig. [31i) the displacement is 
significantly smaller than that of vortices which are part 
of an aggregate (Fig. [3)d) . This confirms that within an 
aggregate the dominant force is the vortex- vortex interac- 
tion (Fint), rather than a force exerted by currents, which 
should affect isolated vortices in the same way as those in 
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FIG. 3. Differential images of typical aggregates of vortices 
before and after a homogeneous heating pulse. White and 
dark spots mark the positions of the vortices before and after 
the pulse. The scale bar for all the panels is the same, (a) an 
aggregate heated to Tmax — O.QSTc. (b) an aggregate heated 
to Tmax — 0.97rc. (c) an aggregate heated to Tmax — 0.99rc 
(d) isolated vortices heated to Tmax ^ 0.97Tc. 



aggregates. Along the same lines, once the vortices have 
dispersed, the intervortex repulsion Fi„i decreases to the 
point where even if we heat the film again the vortices 
do not move. 

When inspecting images taken at T = O.QTTc (Fig.[3lD), 
we can associate individual vortices at their final posi- 
tions with a particular initial aggregate, while ai T — 
0.99Tc this is not possible (Fig. [St). We therefore use 
only the data at T = 0.97Tc for quantitative analysis. 
Even within a specific aggregate, it is not straightfor- 
ward to establish which of the vortices in the initial state 
is which in the final state. We therefore chose to evaluate 
the probability that a vortex is displaced by a distance 
r, assuming that the displacement is minimal. We start 
with a differential image of the whole field of view which 
typically includes ~200 vortices. We then measure all 
the distances between each of the the vortices in the ini- 
tial state and all those in a final state. From this data 
set, we take the shortest distance connecting two vor- 
tices, one in the initial state and one in the final state. 
These two vortices are then removed from the list and 
the process is repeated. This step is done to avoid dou- 
ble counting. The results were averaged over 8 images 
with ~1600 vortices in total. Dividing by the total num- 
ber of vortices gives the probability density Pexp{r)- The 
algorithm works well as long as the the vortices do not 
move further than a typical distance between the initial 
aggregates which is true for data obtained at T = 0.97Tc. 

To compare the experiment with theory, we generate 
images similar to Fig. [3] by numerical integration of the 
equation of motion (Eq. 1). We start with initial config- 
uration of typically 200 vortices taken from experiment 
(e.g. fig. [5]). The integration was performed simultane- 
ously on all the vortices. The time dependence of the 
temperature during the integration was taken from the 



experiment (shown in the insert of Fig. [T]). At each step 
Pint was calculated using the instantaneous positions of 
the vortices. Inside an aggregate, the initial value of Pint 
is of the order of 10~^^N. The pinning force at each tem- 
perature Pp was taken from our transport measurements. 
Pp shown in Fig.l, is consistent with previous measure- 
ments on Nb filmsi^. At T 0.95Tc, for our film thick- 
ness Fp ~ 6 • 10~-'^^iV, more than order of magnitude 
smaller than Pint- 

The explicit form of the drag force rjv is given bj*iS 



Here 77(0) is a sample dependent parameter, and v* is a 
characteristic velocity at which the drag force rjv has a 
maximum. The measured maximal drag force as a func- 
tion of temperature is shown in Fig.[T] Both 77(0) and v* 
were determined for our films from our transport data. 
The values are consistent with prior measurements on 
Nb films^S. Combining measured 77(0) and v* with the 
instantaneous velocity v allows us to determine the drag 
force at each integration step. 

The sudden onset of motion at T* is consistent with the 
non linear form of the viscosity Vortices are acceler- 
ated by Pint- As long as the velocity is below v* the drag 
force increases with v and the motion is highly damped, 
limiting the distance traveled to less than 1 fim. At tem- 
peratures above T*, rjv* becomes smaller than Pint- A 
vortex can be accelerated to w > v* , experiences a re- 
duced drag force, and travel distances significantly larger 
than 1 ^m. 

The accuracy of the simulations is limited by the ex- 
perimental uncertainty in the initial positions of the vor- 
tices which set the exact value and direction of Pint - For 
that reason, the simulations can not reproduce exactly 
the experimental images. We found however that the 
range of the motion of vortices which is important to 
define p*'™ (r) is unaffected by this uncertainty. To sum- 
marize, the pinning force turns our to be small, while 
the maximal viscous force rjv* and Pint are comparable. 
We emphasize that with all the forces known, the vortex 
mass is the only free parameter in the calculation. 

Using the initial positions from the experiment and fi- 
nal positions from the simulations we calculated Psim (r) 
by the same algorithm described above. Similarly, 
Psimif) was averaged over 16 different images. The sim- 
ulated and experimental P{r) are compared in Fig. |4l 
The mass per unit length is given in units of electronic 
mass rue per lattice constant a. In simulations using 
large masses {fj,y > 20?ne/a), the vortices move little 
(r < 2^m) and PsimiT) is large only at small distances. 
For very small masses {fiv < iTie/a), the displacement 
can be large and Psimir) is significant at very large dis- 
tances. Only for masses 5 < fiy < 10, Psim{r) is large for 
2 < r < 5/^771, which is in good agreement with experi- 
ment. This agreement suggests that in the temperature 
window close to Tc, Eq. 1 can successfully describe the 
dynamics of vortices. Furthermore, the results are sen- 
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FIG. 4. Experimental (solid circles) and calculated (bands) 
displacement probability distribution of vortices at T ~ 
0.97Tc. Blue band corresponds to vortex masses of 20 m^ja, 
gray band to 5 rrie/a, and red band to 1 me /a. The high 
value of P{r) at small r is largely due to isolated vortices 
which move less than 1 /im. The error bars and the width of 
the bands represent statistical uncertainty. 

sitive to the presence of an inertial term. Taking into 
account all the experimental uncertainties, our final re- 
sult is 1 < /ii, < 20. 

Most of the literature quotes the vortex mass at T=0. 
The first calculation of the vortex mass by Suhl4 relied 
on the velocity dependence of the free energy of the vor- 
tex core. This gave a mass per unit length /i^ Img/a. 
Additional contributions to the mass so defined were dis- 
cussed by several authors^. Another approach is to define 
the mass through the proportionality between the force 
acting on a vortex and its accelerationjii^ii. The resulting 



mass is of the order of the mass of all the electrons within 
the vortex core. To compare with our results, we need 
the value close to Tc. Suhl used the Ginzburg-Landau 
formahsm^, valid near T^. For Nb at T ~ 0.97Tc Suhl's 
mass is /i^ ^ 10~^me/a. In the approach involving dy- 
namic responsei^, the temperature dependence of the 
mass was calculated by Han^. At T ~ 0.97Tc the mass is 
~ 2 ■ lO'^rrie/a. Our result falls in between these values. 
The only other reported measurement^ comes from the 
response of a vortex array to sound waves. In a large 
crystal, the vortices are tangled and the response of the 
array is collective. This may explain why the reported 
mass is 2 orders of magnitude larger even than the largest 
theoretical prediction. In our experiment we use low den- 
sity vortices in a thin film. At low density, vortices are 
not entangled so the dynamic response is that of a single 
vortex. 

It was suggested that the vortex mass depends on the 
"cleanliness" of the superconductor—. Suhl's prediction 
corresponds to the dirty limit, while the dynamic re- 
sponse prediction to the superclean limit. Our sample is 
intermediate between these limits (mean free path (I k, 20 
nm) roughly equals the correlation length (^o ~ 40 nm)). 
In conclusion, we have demonstrated that the inertial 
mass of a vortex is a meaningful concept which appre- 
ciably influences the dynamics close to Tc- "Quantum" 
vortices may indeed exist in the dirty limit of high Tc 
superconductors. 
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